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Abstract. We prove Central Limit Theorems and Stein-like bounds for 
the asymptotic behaviour of nonlinear functionals of spherical Gaussian 
eigenfunctions. Our investigation combine asymptotic analysis of higher 
order moments for Legendre polynomials and, in addition, recent results 
on Malliavin calculus and Total Variation bounds for Gaussian subor- 
dinated fields. We discuss application to geometric functionals like the 
Defect and invariant statistics, e.g. polyspectra of isotropic spherical 
random fields. Both of these have relevance for applications, especially 
in an astrophysical environment. 
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1. Introduction 

1.1. Motivation. Much effort has been recently devoted to the analysis 
of spherical Gaussian eigenfunctions (to be defined below). These random 
fields are the Fourier components in spectral representation expansions for 
general spherical Gaussian fields, see for instance [I], [9], [13], and [8], [11] for 
extensions; in view of this, their study is of obvious relevance in connection 
with the statistical analysis of spherical data. Namely, the analysis of these 
components and their polynomial transforms (the so-called polyspectra) is 
now one of the leading themes in modern Cosmology, in particular in the 
growing area of Cosmic Microwave Background (CMB) data analysis: we 
refer for instance to [2], [7], [10], [13] and the references therein. 

In short and somewhat vague terms, CMB represents a relic electro- 
magnetic radiation from the so-called age of recombination, e.g. the cos- 
mological epoch when free electrons and protons combined to form sta- 
ble hydrogen atoms; this is now reckoned to have occurred some 3.7 x 10 5 
years after the Big Bang, i.e., some 1.3 x 10 10 years ago. As such, CMB 
radiation is universally recognized as a goldmine of information on pri- 
mordial epochs, and its analysis has drawn enormous theoretical and ex- 
perimental efforts: we refer for instance to http://map.gsfc.nasa.gov/ and 
http://www.rssd.esa.int/Planck for the most relevant ongoing experimental 
activity. 

Among many, one of the leading current themes is the analysis of non- 
Gaussianity of CMB data; for this purpose much effort has been dedicated 
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to the investigation of the properties of the polyspectra mentioned above, 
e.g. polynomial transforms of sample spherical Fourier components. In 
this paper we provide a full characterization for the asymptotic behaviour 
of these transforms, proving Central Limit Theorem results under rather 
broad assumptions. 

From a different perspective, the analysis of high frequency (or high en- 
ergy) Laplace eigenfunctions is a classical subject in Analysis and Mathe- 
matical Physics with much progress in understanding, both rigorously and 
conjecturally, in the recent years. According to Berry's universality con- 
jecture [3], the deterministic high energy eigenfunctions on "generic" mani- 
folds are represented by random monochromatic waves, an equivalent planar 
model to our spherical Gaussian eigenfunctions. 

One in particular interesting aspect of those is their nodal structure [5] 
(e.g. the zero set and the number of its connected components or ovals), 
especially in light of their conjectural relations to percolation and SLE |6j. 
Some results on the geometry of the nodal structure include studying the 
number of the ovals (equivalently the number of nodal components i.e. con- 
nected components of the zero set complement) [16) . and the total length 
of the nodal line [26]. The so-called Defect (i.e. the difference between 
"hot" and "cold" regions for the eigenfunctions, see below or Section [4.2p of 
spherical Gaussian eigenfunctions was addressed in our earlier work |15j . 

1.2. Statement of the main results. Let A§2 denote the usual Laplacian 
on the 2-dimensional unit sphere 5 2 , and consider the Gaussian random 
spherical eigenfunctions fi, I £ Z>o, satisfying 



It is well known that /; can be expanded into the orthonormal basis of 
spherical harmonics {^m} 



(where the coefficients a/ m Gaussian i.i.d.). The random field fi is isotropic 
or rotation invariant, meaning that for any rotation g G SO(3) the dis- 
tribution of /;(•) is equal to the distribution of fi(g ■) (equivalently, for 
every k G N and xi,...,Xk G S 2 , the distribution of the random vector 
(fl(xi), ■ ■ fi{xk)) equals to the distribution of (f t (g • x x ), . . . , f t (g ■ x k ))). 
Equivalently, /; is centred Gaussian, with the covariance function 



where (•,•) is the Euclidean inner product, so that $(x,y) := arccos((a;, y)) 
is the angular distance; here, Pi : [—1,1] — > M. are the usual Legendre poly- 
nomials, see the Appendix for more details. 

The purpose of this paper is to study the asymptotic behaviour of random 
quantities such as 



A §2 f t = -1(1 + 1)/, , E [Mx)} = , E [Mxf] = 1. 




m=—l 



n {x,y) := E[/,(x) • My)] = Pi((x,y)) , x,y G S 2 , 
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where H q {.) are the Hermite polynomials satisfying the differential equations 

H' m (x) = mfl- m _i(x) , E [H m {Z)\ =0,Z~ M(0, 1) . 

As mentioned above, these nonlinear transforms are of interest by them- 
selves as statistical functionals, in connections to the analysis of angular 
polyspectra of spherical random fields (more details to be provided below). 
Importantly, these statistics are the basic building blocks for the analysis 
of general nonlinear functionals of Gaussian eigenfunctions. More precisely, 
it is a well-known general fact (see for instance [19]) that the I? space of 
square-integrable nonlinear transforms of Gaussian eigenfunctions can be 
expanded (in the L? sense) as 

(1) 

G(fi) = f; J -^H q {h) , E [G 2 (fi)} < oo , J q {G) := E • 

As a consequence, the analysis of averaged statistics of the form 

Q{fl) = [ G(fi)dx 
Js 2 

for "generic" G will directly follow from Central Limit Theorem results on 
{hiq\ . To establish the latter, we first need to investigate the asymptotic 
behaviour, as I —> oo, for the variances Var(hi- q ). Note that when both q 
and I are odd, the statistics {hi- q } are identically zero for the symmetry 
properties of Legendre polynomials, e.g. 

Pt(x) = (-l)'PK-*), 

whence integrals of odd polynomials over the sphere are identically zero. To 
simplify the discussion, throughout the sequel we shall consider all limits 
only for even multipoles I. Under this condition, the asymptotic behaviour 
of these variances was investigated in [T5], where it was shown that, for q = 3 
and q > 5 

rn/2 

Var(h q -i) = (4ir) 2 q\ if(costf) smtfdtf ~ (47r) 2 g!-f , 
Jo ' ' 

/•oo 

c q = / iPJoW q diP > . 

J 

For q = 2,4, the order of magnitude of the corresponding variances is larger 
(see Lemma [33]) : 

Var(h q ,) « { 4/ ^ " = 2 4 . 

The constants c q are immediately seen to be strictly positive for all even 
values of q. For odd values, we conjecture this to be always the case 5 a 
formal proof is left for future research. Therefore, the statement of our 
main result will entail this condition explicitly: 



Theorem 1.1. For all q such that c q > 0, we have 

/lr>;.„ 
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Theorem 11.11 is a building block for a more general claim: under minimal 
regularity conditions (i.e. the existence of at least one nonzero coefficient 
J q {G) corresponding to a nonzero term c q > 0) we shall have a CLT for 
square integrable nonlinear functionals of spherical Gaussian eigenfunctions, 
i.e. for any Q(fi) = f g2 G(fi)dx such that EG 2 (fi) =: <Tq{1) < oo, 

g(/2i)-E[g(/ 2t )] 

-nzr, >dJv{0, 1), as I ->■ oo. 

In fact we will prove the following stronger result: for q ^ 4, 

(2) ri TV f /T/ ^ . ,^(0,1)) =0(l^ 

for some 6 q > 0, where c£tv(-j .) denotes the usual total variation distance 
of random variables 

d TV (X,Y) = sup \Pv(X €A)- Pr(Y e A)\ , 
AsB(R) 

with the Borel a- field B(R). For q = 4 the rate of convergence ([2]) is of 
slower logarithmic order. 

1.3. On the proofs of the main results and outline of the paper. 

The ideas behind our main argument can be summarized as follows. Be- 
cause fi(-) is a Gaussian field, for any fixed x G § 2 , H q (fi(x)) belongs to 
the so-called q-th order Wiener chaos generated by the Gaussian measure 
governing /$(•) (see [E]), and so does any linear transform, including hi ;q . 
As a consequence of the Nourdin-Peccati Theorem for Stein-Malliavin ap- 
proximations of Gaussian subordinated random variables (see for instance 
Theorem 5.26), the following bound holds for each even I, q > 2 : 



(3) d TV - h M f ,AA(0,1)| <2 l q - 1 ( cum ^n) 



K y/Var(hi. g Y " ' ') ~ V 3g \Var 2 (h l;q ) 

here, cum^ (Y) is the 4th order cumulant of Y, see [18], [19] for more discus- 
sion on these points. 

The latter bound shows that if we prove that 

(4) cum A (hi- q ) = o^oo (Var 2 (hi- q )) , 

the Central Limit Theorem for hi- q (where q is fixed and I — > oo) will follow. 
The bound (jl]) for q > 5 is proved along Section [2j Here we first express the 
4th order cumulant as an integral over (5 2 ) 4 , using the well-known Diagram 
Formula (see Section 12. 1|) . This will allow us to obtain the desired bound (jj]) 
via a tricky multiple application of the Cauchy-Schwartz inequality (Propo- 
sition 12.21 whose proof in Section 12.31 takes on most of Section [2]) ; to this 
end we will divide the domain of integration into the "local" and "global" 
ones (for the definitions of various ranges see Section [232]). A more detailed 
explanation of the proof of Proposition 12.21 may be found in Section 12.3.11 

For q = 3, 4 proving the bound (jU) will require special arguments, pre- 
sented in Section [3l While the case q = 3 was already covered earlier in 
[12] . for q = 4 the asymptotic analysis requires the evaluation of so-called 
Gaunt integrals, connecting moments of Legendre polynomials to Wigner's 
and Clebsch-Gordan coefficients (see [13] 
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Various applications of Theorem 11.11 for the Defect statistics and general 
polyspectra are discussed in Section [H The basic idea is to use the Central 
Limit Theorem for h[. q to establish asymptotic Gaussianity of finite-order 
generic polynomial sequences, and then exploit the expansion (TjQ) of arbitrary 
functional (e.g. the Defect) to prove that the distribution of any such 
functional can be asymptotically approximated by means of a finite-order 
expansion. 

1.3.1. Some conventions. In this manuscript, given any two positive se- 
quences a n ,b n we shall write a n « b n if there exist two positive constants 
ci,C2 such that c\a n < b n < C2d n , for all n = 1,2,..., and a n ~ b n if 
lirUji-^oo a n /b n = 1. Also, we shall write a n <C b n or a n = 0(b n ) when the 
sequence a n /b n is asymptotically bounded; we write /J,(dx) for the usual 
Lebesgue measure on the unit sphere, so that J* g2 fJ,(dx) = Air. 

1.4. Acknowledgements. The authors would like to thank Giovanni Pec- 
cati, Zeev Rudnick and Mikhail Sodin for many stimulating and fruitful 
discussions about the subjects raised in the present manuscript. A sub- 
stantial part of this research was done during the second author's visit to 
University of Rome "Tor Vergata", and he would like to acknowledge the 
exemplary hospitality of the institution, and the generous financial support. 



2.1. Some preliminaries. We shall now focus on fourth order cumulants 
for Hermite transforms for arbitrary q > 5. First we need to recall some 
well-known background material on the so-called Diagram Formula (see for 
instance [21] . |19j . |13| or |17] for recent textbook references). 

Fix a set of integers a±, ...,a p . A diagram is a graph with {a\ + ... + a p ) 
vertexes labelled by l,...,p, (ot\ vertexes are labelled by 1, 02 vertexes are 
labelled by 2...) such that each vertex has degree 1, i.e. the edges have no 
common endpoints. We can view the vertexes as belonging to p different 
rows and the edges may connect only vertexes with different labels, i.e. there 
are no flat edges on the same row. The set of all such graphs 7 is denoted 
by r(«i, a p ); we write Tc{ot\, a p ) for graphs that are connected, i.e. 
it is not possible to partition the vertexes into two subsets A and B such 
that no edge connects a vertex in A with one in B. 

Given a diagram 7, let 



i,3 

and, moreover, the multiplicities of opposite edges equal, as the following 
lemma states. 

Lemma 2.1. Let 7 G Tc(q,q,q,q) with arbitrary q > 1, and r\ = 77(7). Let 
e = (i,j) any edge in 7 and e' = the unique edge with vertexes disjoint 

with e, so that {i,j,i',f} = {1,2,3,4}. Then rj e = r) e i . 



2. The central limit theorem for hi ;q , q > 5 



2(7) = fai(7)) e z(S) 
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Proof. The statement of the lemma follows immediately from the fact that, 
by reordering the vertexes corresponding to the same label (i.e. in the same 
row) if necessary, we may assume that all the edges are between vertexes in 
the same column. □ 



2.2. Cumulants of hi- q . With notation as above, the well-known Diagram 
Formula (see e.g. (19j ) provides the following neat expression for the cumu- 
lants of our statistics {hi q } , namely 



= J cum{H qi (fi{xi)),...,H qp (fi(x p ))} dxi...dx p 

(S 2 )p 

= j ^2 Y\ E iM x i)fi( x j)} dxi...dx p 

e / n(v{fi(x i )f l (x j )}f^dx 1 ...dx p , 

7erc(giv,<?p)( S 2)p 

(recall that 77(7) is the vector of multiplicities of edges in 7 as above). 

We constraint ourselves to the 4th order cumulants of the form cum(hi ;q , hi- q , hi- q , hi- q ); 
in this case the expression above simplifies to 

(6) cum{hi. q ,hi. q ,hi. q ,hi-,q) = E M (v(l)), 

7er c (<?,<?,<?,?) 

where for a vector r\ G Z> we set 




with the shortcut dx = dx\ ■ . . . ■ dx^. Since © is a finite summation of the 
C{rj) with rj corresponding to a connected diagram 7 £ Fc(q,q,q,q) it then 
remains to bound each of terms as above separately. 
We know [13], Lemma 5.2, that for q > 5, 

Var(hi. q ) ~ y| 

and aim at proving the bound 

cum(ht q , hi. q , ht q , hi- q ) = o 

(or stronger); this is sufficient for the central limit theorem by O The 
following proposition implies this estimate bearing in mind ([6]). 




ON NONLINEAR FUNCTIONALS OF RANDOM SPHERICAL EIGENFUNCTIONS 7 



Proposition 2.2. For q > 5 and 7] = 77(7) with 7 a connected diagram, one 
has 



\M(rj)\ = 

In particular, for every I > 5, 



O(^f) q = 6 



\M(r L )\=o(j £ 

The methods of the present section also give a useful (i.e. smaller than the 
square of the variance, sufficient for central limit theorem) upper bound for 
M(rj) with r] corresponding to q = 4. It is however weaker than the precise 
asymptotics of Lemma 13.31 which is the reason for a special treatment we 
gave to q = 4. 

Remark 2.3. In fact, it has been recently established by G.Peccati and 
coauthors that only a subset of terms as above need to be bounded to es- 
tablish the CLT - those corresponding to circular diagrams (i.e. diagrams, 
all of whose rows are linked with precisely two other rows, see e.g. |19| . 
Proposition 11.2); the latter were easier to bound by our earlier methods. 
We will however not use this observation as our present methods can cope 
with arbitrary terms, though resulting in slightly weaker upper bounds. 

2.3. Proof of Proposition [2T2l 

2.3.1. On the proof. One observes that for a "generic" point, each of the 6 
terms Pi((xi,Xj)) in the integrand product in ([7]) is bounded by -K as in 
(|33|) , Therefore, unless at least one of the 6 angles involved is small, the 
integrand is of order O fw] , and the total contribution to the integral in ([7]) 
for q > 5 should be of order smaller than 75, sufficient for the CLT0. 

To quantify the latter statement, we introduce a small parameter e = e(l) 
and separate the domain of integration (S 2 ) 4 into the set C(e) of points 
x G (5 2 ) 4 all of whose angles ${xi,Xj) are greater than e, and its comple- 
ment ($(xi,Xj) < e for least one pair of indexes (i,j)), see Section T2.3.2I 
We call the contribution of the latter subdomain "local" and, analogously, 
the former's contribution, "global". The global and local contribution are 
bounded in sections 12.3.31 and 12.3.41 whence it remains to choose the optimal 
parameter e the arises from the tradeoff we get from those bounds (Section 

To bound both the global and local case one employs the following ob- 
servation. It is possible to decrease the number of different angles $(xi,Xj) 
involved in the integral 12.3.21 by applying the Cauchy- Schwartz inequality, 
multiple times if necessary. The upshot is that in integrals like 



/ 



P z (cos e((x 1 ,x 2 )) Sl Pi(cos e{{xz,x A )) S2 dx 



^This heuristics is not entirely correct, as we believe the correct order of magnitude of 
the 4th order cumulant to be proportional to rather than ^ . It means that the regime 
where at least one angle is small does contribute to the integral. 
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the variables split and we end up evaluating moments of individual Legendre 
polynomials (that were readily evaluated, see (fT2j) and (fTT|) ). Therefore, 
when we apply the Cauchy- Schwartz inequality to reduce the number of 
different angles, it is beneficial to pair up angles corresponding to disjoint 
edges in the diagram. For q>7 these kind of observations, combined with 
uniform estimate (|14p valid on C(e) and the small measure of its complement, 
are sufficient for our purposes. The cases q = 6, 5 are a bit more subtle as 
in these cases we will have to exploit the special structure of vectors n 
corresponding to connected diagrams. 

2.3.2. Global and local terms. Recall that for a rj £ Z> we defined C(rj) as 
in 0. For a small parameter e = e(l) > we decompose the domain of 
integration in ([7J as following: 

(S 2 ) 4 = C(e)UC(e) c , 

where 

C(e) := {(xi,... x 4 ) G (S 2 ) 4 : Vi < j, 0(x f , Xj ) > e} 
and C(e) c is its complement. We then write 

M {rf) = M glob {rr, e) + M [oc {ry, e), 

where 

(8) M glob (ms)= f\{P{{x u x 3 ))^dx 

C *<3 

is the "global contribution" and 

M loc (ry,e) = J J] P({x it Xj ))^dx 

is the "local contribution" . 

The global and local contributions are bounded in the following couple of 
lemmas, whose proofs are given in Sections 12.3.31 and 12.3.41 respectively. 



Lemma 2.4 ("Global contribution"). As I — > oo, q > 5 arbitrary, and 
V = ^(7) with 7 G Tc(q, q, q, q), we have 

q = 5 

\M glob {ry,e)\ = lO{^f) q = 6 




liei~ 3 



q>7 



with constant involved in the "O" -notation depending on q only. 

Lemma 2.5 ("Local contribution"). Under the assumptions of Lemma \2.J\ 

and the additional assumption 

(9) e » j, 

one has 

'(log/) 3 /VV 2 \ 



(10) \M loc (rj;e)\=0 



P 
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To prove lemmas 12.41 and 12.51 we will use asymptotics for the 2nd and 4th 
moments of Legendre polynomials: as I — > oo 



(11) / Pl(t) 



l 

2 dt ~ C 2 y 





for some C2 > (cf. ([5T]) ). and 



i 



(12) J m) 



4r. fog I 

at ~ C4- 



P 





for some C4 > (cf. Lemma l3.2|) . 
2.3.3. Bounding the global contribution. 



Proof of Lemma \2.4\ By the definition of the domain of integration £(e), 
all the angles involved satisfy "d(xi,Xj) > e. Let 77 = 77(7) corresponding to 
a connected diagram 7 G Tc(q, q, q, q). 

First we assume that q > 7. By an easy counting argument, ([5]), and 
Lemma [2 .H it may be shown that 77^ > 2 for (at least) two disjoint elements 
of 77. With no loss of generality we assume that both 7712 > 2 and also 7734 > 2. 
Since 7 is connected, we may also assume 7713 > 1, 7/24 > 1, again, by the 
virtue of Lemma 12. 11 

Consider the integrand 

(13) K ri {x u ...,x±) = Y[P l ({x i ,x j })^ 

i<j 

of ([8|). On C(e) for every K j we have the uniform upper bound 

1 



(14) |P / ((x,,x J ))|« 
by (USD- Hence 

\Krjix!,.. . ,x 4 )| < ^ £ | g _ 3 i 3 K(^i^2)) 2 ^K(^3,^4)) 2 |^((a:i,3:3))-PK( a; 2,^4))|, 

so that 
(15) 

|M gW (77_;£)| < / J Pz((xi,X2)) 2 P/((x 3 ,x 4 )) 2 |PK^i^3)) J P/((2;2,2;4))Mx 



(Ze)?- 3 



(note that at this stage we can afford to increase the domain of integration 
to the whole of (S 2 ) 4 ). 

In order to treat the latter integral we apply the Cauchy- Schwartz in- 
equality, dividing the 4 terms into pairs. We team up each edge with its 
disjoint complement (i.e. the unique edge with no common vertex); it is 
then possible to split the variables in each of the resulting integrals. This 
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approach yields: 

Pl((x 1 ,X 2 )) 2 P l {(x 3 ,X A )) 2 \Pl((x 1 ,X 3 ))Pl((x 2 ,X 4 ))\dx 
(§ 2 ) 4 

/ \ 1/2 / x 1/2 



< 



Pi({xi,x 2 )) Pi((x3,X4}) dx 



\(S 2 ) 4 
log; 



P z ((xi,x 3 )) 2 Pz((x 2 ,x 4 )) 2 dx 



by (|12p and (|lip . The statement of the present lemma for q > 6 then follows 
upon substituting the latter bound into (fl"5j) . 

For g = 5,6 the same argument remains valid; however the bound it gives 
is insufficient, and we will need to exploit the special structure of r] in these 
cases. For q = 5, 77 has one of the following three shapes (in some order): 

(16) 7? = (4,4,1,1,0,0) 
or 

(17) r? = (3, 3, 2, 2, 0,0) 
or 

(18) r? = (3, 3, 1,1, 1,1), 

where, by Lemma [27TT the edges corresponding to rjij = 4, 3, 2 are disjoint. 

For the first case (|16p . we have (up to reordering {xi}) 



Mgi ob {rr,e)= J Pi({x 1 ,x 2 )rPi({x 3 ,x 4 )rPi({x 1 ,x 3 ))P l ((x2,x^)dx, 

C(e) 

so that the uniform bound (|14p yields 

\M glob (rr,e)\ < i j J P l ((x 1 ,x 2 ))*Pi((x 3 ,x 4 ))Ux^^^ 

C(e) (S2)4 

by (fT2|) . Next, for the second case (fT7|l . we bound 

\M glob (rr,e)\ <-jiJ Pi^xux^fPi^x^x^fPi^x^fP^x^x^fdx, 
(s 2 ) 4 

and continue as for q > 7 by bounding the latter integral using Cauchy- 
Schwartz, pairing together disjoint edges. 

Finally, for r\ as in (fl~8|) . we use (fT4"P to bound 



< r / ^K(^1^2)) 2 i 3 K(^3,^4))VK(^1^3)) J PK(^2,2;4))PK(^1^4))fl((x2,2;3))Mx. 



|^o&(f?; e)| 

1 

(§ 2 ) 4 

Pairing up the 4th degree terms in the latter integral together and the rest 
separately, we apply Cauchy-Schwartz as before, twice for the second term. 
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This leads to 



Pl({x 1 ,X2)) 2 Pl((x 3 ,X^)) 2 \Pl({xi,X 3 ))Pl((x2,X 4 ))Pl({xi,X4))Pl((x2,X 3 ))\dx 



(S 2 ) 4 



< 



W) 



P l ((x 1 ,x 2 )) 4 P l ({x 3 ,x i )) 4 dx 



Pi({xi,x i )) A Pi((x2,x 3 )) A dx 



1/2 



1/4 



mx^x^fp^x^fdx 



W) 



this implies the statement of the present lemma for q = 5 via separation of 
variables and (|12p . The proof for q = 6 is very similar to the above (but 
somewhat easier) and thereupon omitted here. In this case r\ is of one the 
following 5 forms: 

7? = (5, 5, 1, 1, 0, 0), (4, 4, 2, 2, 0, 0), (4, 4, 1, 1, 1, 1), (3, 3, 3, 3, 0, 0), (3, 3, 2, 2, 1, 1) 

□ 

2.3.4. Bounding the local contribution. 

Proof of Lemma \2.h\ We may assume with no loss of generality that 

■&(xi, x 2 ) < e 

in the relevant domain C(e) c . We divide the possibilities for r\ into three 
cases, up to reordering (that r] falls into one of those cases follows from © 
q > 5, Lemma I2TT1 and the the connectedness of 7): 

(1) For every i < j, rjij > 0. 

(2) We have 



V12 = ??34 > 4, 7/13 = V24 > 1, f?14 = ??23 = 0. 



(3) We have 



?7l2 = ??34 > 2, r/13 = 7/24 > 2. 



In case Q] (which gives rise to the dominating term) we may use Cauchy- 
Schwartz twice to bound 



i<j 



(19) 



\M loc (rr,e)\< J W\P{{xi,Xj))\dx 
( 



1/2 



< 



^((^1^2) 2 i 3 K(^3,^4) 2 -P/((a ; 1^3) 2 ^ X 



1/2 



P\ ( (xi ,X4 ) 2 Pi ( (x 2 , x 3 ) 2 Pi ( (x 2 , x± ) 2 dx 



V c ( £ ) 
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and bound each of the two integrals of (|19|) separately. For the first integral 
we use Cauchy-Schwartz again: 



Pi ( {x 1 , x 2 ) 2 Pi{ (x 3 , Xi ) 2 Pi ( (xi , x 3 ) 2 dx 



1/2 



1/2 



< 



Pi((x\, x 2 ) 4 dxidx 2 Pi((x3, x^dx^dxi 



W) 

log/ (log/) 1 / 2 (log/) 3 / 2 



\{x2:i?(xi,a;2)<£} 



Z 2 I 



by separation of variables and 



({x 2 : x 2 ) < e}) < 



For the other integral of (|19p . we may use the lack of symmetry w.r.t. 
variables to improve the bound as follows: 



Pi((xi, x 4 ) 2 P[({x 2 , x 3 ) 2 Pi({x 2 , Xifdx 



1/2 



< 



Pi((xi, x^) 4 dxidx4 



1/2 

Pi((x 2 ,x 3 ) 4 dx 2 dx 3 I x 



\{x 2 :i3(; 

(log/) 3 / 2 



r V 

;i,a'2)<e} / 



\{a;2:'!?(^i i ^2)<e} 
1/2 



(log/) 1/2 (log/) 1/2 (log/) 1 / 2 
Z ' I £ ' I £ 



e 2 . 



We then obtain the statement of the present lemma in this case upon sub- 
stituting the last couple of estimates into (|19|) . 
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In case [2] we use similar ideas (Cauchy- Schwartz twice) to obtain 
M loc { m e)\< J P l ((x 1 ,x 2 )) 4 P l ((x 3 ,x 4 )) 4 P l ((x u x 3 ))P l ((x 2 ,x 4 ))dx 

( \ 3 / 4 

Pl{(x 1 ,X 2 )) 4 Pl((x 3 ,X4)) 4: Pl((x 1 ,X 3 })Pl{(x2,X4})dx I X 



< 



W) 

( 



\ 



1/4 



P 1 ((X 1 ,X 2 )) 4 P 1 ({X 3 , X 4 » 4 ^«X1, X 3 )) 4 ^«X2, X 4 » 4 dx 



(log/) 3 / 2 

/ 3 

(log/) 2 



/ 



1/4 



P^zi,^)) 4 ^^ 



\(§ 2 ) 



Pi({x 2 , X4)) 4 dx 2 dX4 : 



\{x2--'&(xi,X2)<e} 



/ 4 



which is smaller than the RHS of (|10p by ([9]). 
Finally, for case [3] we similarly have 

(20) 

\M l0C (rr,e)\ < J \P l ({x u x 2 )f\\Pi({x 3 ,X4)) 3 \Pi((x 1 ,x 3 )) 2 Pi((x 2 ,X4)) 2 dx 
< J mx^x^fP^^x^fmx^x^fPiiix^x^fdx 

me) 

1/2 / x 1/2 



< 



Pi((xi,x 2 )) i dx 1 dx 2 



us 2 ) 2 



Pi{{x 3 ,Xi)) A dx 3 dx4 



1/2 



Pl{(x 1 ,x 3 )) A dx 1 dx 3 



\ 



1/2 



Pi((%2,za)) dx 2 dx4 



\(S 2 ) 2 / \{x 2 :f(^i^2)<£} 

<<; log/ (log/) 1 / 2 (log/) 1 / 2 ^ = (log/) 2 c 



/ 



Z 2 



/ 



/ 



/ 4 



which is less than latter of the expressions on the RHS of (fTU|) . again by 
®. □ 



2.3.5. Concluding the proof of Proposition 001 

Proof. In order to finish the proof of the present proposition it remains to 
make a suitable choice of the parameter e(Z) so that both the local and the 
global contributions as bounded by Lemmas 12. 41 and 12. 51 will be smaller than 
the expressions prescribed in Proposition 12.21 The optimal choice for the 
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arising trade-off is 



s(l) 



1 
i 

t G/7 



q = 5 
q = 6 
q>7 



giving the bound in the statement of the present proposition. 



□ 



3. The Central Limit Theorem for hi- q , q = 3, 4. 

We shall start from the investigation of total variation bounds for hi-s; 
this result was established in [T3], see also [T2], [TJ], but nevertheless we 
report it here for the sake of completeness. The Lemmas below make some 
use of so-called Wigner's and Clebsch-Gordan coefficients; see the Appendix 
for their definition and discussion of some important properties. 

Lemma 3.1. 1. The variance of hi-3 is given by 

(21) E K;3 ]=6X(4^(< < <) 2 ~^£ ; 

2. For the fourth-order cumulant of hi ; s we have 

(22) cum 4 (hs) - -IpY ( o 

3. The following total variation bound holds: 

As argued in the Introduction, the general strategy for the proofs of our 
convergence results requires a careful evaluation of the variance and suit- 
able bounds on fourth-order cumulants. For q = 4, these computations are 
provided in the two Lemmas to follow. 

Lemma 3.2. (1) The variance of hi-4 is given by 

i 

Var{h tA } =4\(4ir) 2 J P^dt 
o 

(2) As I — > oo we have 

1 mfdt 3 lORl 



2tt 2 I 2 ' 







In particular, 



Var{h t . A } ~ 24 2l ° gl 



I 2 

Proof. Since 

E[/l,; 4 ] = 0, 
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the first part of the present lemma follows from 



H 4 (fi(x))dx ■ / H A {fi(y))dy 



Var{h l . A }=E[hf. 4 \=E 

Ls 2 

l 

E [# 4 (/^))#4C/Ky))] dxdy = 4! J P^x, y)fdxdy = 4!(4vr) 2 J P,(/ ) ',// 

§ 2 xS 2 S 2 x§ 2 

To see the second part, we invoke Hilb's asymptotics f)32[) in the Appendix 
to write (up to an admissible error, as it is easy to directly check) 

1 1+1/2 I 

J Wdt ~ ~ J sin (jq^) ^oW) 4 # ~ J / J • W) 4 # 

o l 1 

1 /" 4 sin(^ + vr/4) 4 1 3 f dip 



/ 5 



i i 

by the standard asymptotics for the Bessel Jo- □ 
Lemma 3.3. As I —> oo, we have 

cum^ {hi-^} ~ Z -4 . 



Proof. For our purposes, we need to show that 
(24) 

A 1 := / P l ((w,z))Pi\(w,w'))Pi((w',z'))Pj i ((z',z))dwdzdw'dz' = O ( 1( ' K " ' 

Js 2 x....S 2 \ 

and 
(25) 

A 2 := / Pf(( W ,^))Pf« U ;, U ; / »P J 2 « U ; / ,/»Pf«^,2;»d U ;dW^ / « ^ 

Js 2 x....S 2 1 

Concerning the first term, we note that 

/ Pi{{w,z))P?{(w,w'))dw 
J§ 2 

f A \ 4 



x / Yi mi (w)Y lml (w)Y lm ,(w)Y lml (w')dw 



(21 + 1) 2 S y «mi( 2; ) y «m^( w/ ) y K(w / ) y /m^(^ / ) 



3 
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qLM (jL,-M 

x l^^ 1 ) i°io/o) 2L + 1 

LM 



Likewise 



(2 ( ^\ 3 )2 E Y lm ,(w')Y lm2 (z')Y lm4 (z')Y lm4 (z' 

iL'M' r>L',-M' 



L'M' 



E(-!) L+M C ' 



^0/ 2L/ + 1 



Note that L is necessarily even here, otherwise the Clebsch-Gordan coef- 
ficients {C{$ } are identically null from flM|), whence (-l) L+M = (-1) M . 
Iterating the same argument twice more we obtain 

(jLM qL-M qVM 1 qL',-M' 

2^ zJ _1 J i G «o/oi Z^ _1 J \ G «o/oJ 2L' + 1 

mi...mJLM L'M' 

(26) 

(jL"M" sjL",—M" fjL'" M" sjL'" ,—M" 

x \°w/o/ 2L" + 1 ^ ^ ' \°«o/o / 2Z7" + 1 

L"M" " " L"'M"' 

Now, applying iteratively the orthogonality identity (|37p . we have 

(jLM (jL-M qL'M' qL',-M' 

Z_> l^\- l > 2LTI 1^020/ 2L> + 1 

m' 2 rn' 3 m' 4 IM VM> " ! 



\ -> \ -> X^MOJ Irrnlm'z m\m' 2 _ \ - {^MOS x m 'l 

2-f ff f 2L + 1 2L + 1 " V ( 2/ + !)( 2L + 1) mi ' 
Applying the same argument to the last two terms in (|26|) . we obtain that 

MS,} 4 {c»} 4 C! 



2L + 1 2L' + 1 (2Z + l) 2 

l T.T.I y ' 



J; / {^o°o} 1 

+ 1 1 r 2L+1 / 



4 ^ 2 , 9 

1 1 log 2 / 



o 



21 + 1 I ^ 2L + 1 [ V 2Z + 1 Z 4 

Here it is interesting to recall that 

see also |14j . Lemma Al and the proof of Lemma 2.3 therein. 

Let us now focus on (J25j) . Using again (p£2|) and (jHJ), we have that 



P 2 (( U ;,z))P 2 (( U ;, U ; / »^ 

2 



ON NONLINEAR FUNCTIONALS OF RANDOM SPHERICAL EIGENFUNCTIONS 17 



4?r 



2Z + 1 J J § 2 



f 4 1 4 

1 2/TT | ^ Ylmi ^ Ylm ' 2 ( Z "> Ylm ' S ( W ') Y lm' 4 (w) 

mim,2m' s m' 4 

x / Y lmi {w)Yi m2 (w)Y lm > ' (w)Y lm > (w')dw 

(4-^)3 ^ ^ 

^ + jj2 £ y im 1 (^)^/m 2 (2)^ m ; ! ('W / )^m^(^ / ) 



2L + 1 

LM 



Iterating the argument, we find that 

m - E E(-!) M {«} 2 ^gf^ e.(-!) m ' {^;°o} 

(27) 



2L + 1 ^ ' " l~ (Utu J 2L + 1 

L'M' 



sjL"M" qL",-M" pL"M"' n L ,-M' 1 



E (- 1 ) M " {^o'/o } -f" E (- 1 ) A/ "' {«} 

L"M" L"'M"' 

Again it is sufficient to apply (f3T|) four times to have 

E E {C/o?o} {Oozo} 

L,L'L",L"' M,M'M",M"' 



( r L"0\ 2 ( r L"'0\ 



9 zL' sL" sL'" %L sM' xM" sM'" rM 

°L °L' °L" °L'" °M V °M" °M'" 

(2L + 1)(2L' + 1) (2L" + 1)(2L /// + 1) 

8 / , , r- \ 8 



=E(JiS)=D^)(iiS)' 

LM v ' L v y 

We shall now need the following result (see again [II], Lemma A.l): 



I I L\ 2 2 1 



OJ 7iL X vr X L(2/-L) 1 /2(2/ + L)V2 ' 2 - 7iL - 5 
Simple manipulations then yield 

e<-W^)Ve (2L+1) 



J - L*(2l- L) 2 (2l + L) 2 

L — 
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' (2L + 1) ^ (2L + 1) 



L=0 



L 4 (2Z-L) 2 (2Z + L) 2 j^ i L\2l-L) 2 {2l + Lf 



J_ (2L + 1) 1_ 2 ^ (2L + 1 

- A/4 2^ r,4 + 74 2^ 



4/ 4 ^ L 4 Z 4 £-< L(2/-L) 2 (2/ + L) 



2«-2 



L= 



which completes the proof of the upper bound. To prove that this bound is 
sharp, it suffices to notice that 

□ 

Combining the variance and cumulant results, and exploiting ([3]), one 
finally obtains the following result. 

Proposition 3.4. As I — > oo, we have 

d TV [ - hi] \ N ,A/-(0,1)| =0 



yTar(/t 4; /J y \ lo g z 



4. Applications 

4.1. Polyspectra for spherical random fields. Let T{x) be a zero-mean 
Gaussian and isotropic spherical random field, i.e. a measurable application 

T : S 2 x Q, — > M such that = T(gx) for all elements of the group 

of rotations g G SO (3). It is well-known that the following mean-square 
representation holds, in the L 2 (dx x dP) sense (see [13j, Chapter 5): 

T(x) =Y^Ti{x) , where A^T, = + 1)T, . 
i 

We can hence view the eigenfunctions fi as the normalized Fourier compo- 
nents of such spherical field, e.g. fi(x) := T/(x)/-^/e [T 2 (x)]. In this subsec- 
tion, we shall consider the central limit theorem for polynomial functionals 
of the form 

Q , 

Zi = J2b q {fi(x)} q dx , for some Q G N , b q E R . 

9=0 *' 52 

When we view the eigenfunctions /; as the Fourier components of an isotropic 
spherical random field, these polynomial statistics cover, for instance, the 
well-known (moment and cumulant) polyspectra of the random field. These 
are the crucial statistics when searching, for instance, for possible non- 
Gaussian behaviour in T{x)\ see for instance [2], [10] : and the references 
therein. Note that there exist deterministic coefficients /3q,...,/3 p such that 
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we can write 

Q , Q 
z i = /~2P q / H q {f 2 i{x))dx = ^2(3 q h 2 i ]q . 

q=0 g=0 

From the results in the previous Section, we have immediately the following 

Corollary 4.1. Assume that c q > for at least one q such that (5 q ^ 0. 
Then 

Zi-E \Zr] . „ . 

' = AT(0, 1) , as I -> oo . 

The proof is immediate in light of Theorem 11.11 Indeed, we are dealing 
here with a finite linear combination of asymptotically Gaussian random 
variables, and we recall that for random vectors with components in Wiener 
chaoses, the multivariate Central Limit Theorem follows from convergence 
in distribution of the univariate components, see [20J. This result thus ex- 
tends the Central Limit Theorem provided in [12] for the sequence {/i2Z;3} 
to polyspectra of arbitrary orders. 

It is actually possible to establish stronger results, i.e. to study the rates 
of convergence in the total variation bound. Rather than focusing on this 
issue, we move to the Central Limit Theorem for the case of a more general, 
infinite-order 1? expansion, as it is the case for the Defect. 

4.2. Defect. In this subsection, we shall focus on one of the most important 
geometric functionals, namely the Defect. The Defect (or "signed area", see 
|5j) of a function i/j : 8 2 — > R is defined as 



(28) V(ip) := meas (^ _1 (0,oo)) - meas (^ _1 (-oo,0)) = / H(ip(x))dx. 
Here Ti(t) is such that 



(29) W(t) = l[o, o)(*)-1(- O o,o](*) 



1 t > 
-1 t < 
4 = 



where Ia(^) is the usual indicator function of the set A, and dx is the 
Lebesgue measure. In our case, the Defect is the difference between the 
£1X6 clS of positive and negative inverse image of fi, denoted 

T>i :=©(/i). 

It has been shown by [14] that the following expansion holds, in the 
L 2 (dP) sense 

P y h A(-l) (2g-l)ll 

Vl ~ \ (27TT)! hl > 2q+l - ^ V^WTW l ' 2q+1 " 

q=l 5=1 

Trivially E [T>i] = 0, and from [15] we have that 

C2g+1 . n - 2 , ( 2 <?) ! 

q=l 



oo 



Var(V,) = E [Vf] ~ £ a,^ + o(r>) , a, = 
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are the (suitably normalized) Taylor coefficients of arcsin are asymptotic to 



1 



q ~ 2^q 3 / 2 
by Stirling's formula, and 



+ o(q 



-3/2^ 



as q — > oo 



32 



9=1 



27 



Note that we know C3 > from (|21|) : any term corresponding to C2 9 +i = 
can simply be dropped from the expansion, so the rate for this variance is 
precise. In view of Theorem ll.il it is then not difficult to prove the following 
result. 



Corollary 4.2. As I — > 00, we have 

V21 



y/Var{V 2l ) 



^dAA(0,l). 



Proof. The proof follows a standard argument for nonlinear transforms of 
Gaussian measures, see for instance [21]. Define 



<^2q+l 

{ (2? + l)! 



h,: 



using the trivial inequality E [(A - B) 2 ] < 2E [(A - C) 2 ] + 2E [(C - B) 2 ] , 
we have that 



E 



< 2E 



< 



V, 



v,. 



i.rn 



( T^l T^hra 



Var(Vi 



y/Var{V{) y/VarfPt) 



-E [(D, - P /;m ) 2 ] + 2 



+ 2 E 

Far(£>; ;m ) 
Far (D,) 



V, 



l:m 



y/Var(Vi) y/Var(Vi ir 



+ 1 



2y/Var(V, 



I'm ) 



Now, using the same argument as in the proof of Proposition 4.2 from 
[TB] . pages 9-10, we have that 

E^-V,^) 2 ] = E{-^t} 2 IE^2 9+1 ] 

q=rn ^ ^ ' ^ 



1 V 



q=m 



< 



1 1 v-^ C5 , 2 - / 1 

v q=m * \ v 



so that 



Far (A) 



E [(D, - P i;m ) 2 ] = O 
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and 



Var(Vi. m ) i 2y/Var(V l;m ) 
Var{Vi) y/Var(Vi) 



It follows immediately that 



2 + 0{ 



O 



2W1 + 



m 



(30) 



E 



V. 



l:m 



y/Var{Vi) y/Var(V l;m ) 
Now, for every fixed m we have 



O 



1>l;r 



-^■ d jV(0, 1) as i — > oo, 



and since m can be chosen arbitrarily large, the random variables ^ , 

[ \JVariPi) 

must have the same limit, bearing in mind (|30p (see e.g. |21|). □ 



5. Appendix 

The Legendre polynomials are defined by Rodrigues' formula 

Legendre polynomials are orthogonal with respect to the constant weight 
ui(t) = 1 on [—1, 1], indeed 



(31) 



1 25 l2 
P h {t)P h {t)dt=y-^ l - 



they also satisfy the well-known Hilb's asymptotics (see e.g. [23J, formula 
(8.21.17) on page 197): 



(32) 



Pz(cosi 



sin ( 



1/2 



J ((/ + 1/2)6) + 5(9) , 



where Jo is the standard Bessel function, and the error term satisfies 

m « { e T 0{rm) • d " <0< T 

V ' I 9 2 , < 6 < cl-\ 



In particular, for 9 £ [0 
(33) 



' 2J ' 



Pi(cos9) < 



19 



Let us now review briefly some notation on Wigner's 3j coefficients; see 
[25], [23] and [1] for a much more detailed discussion, in particular concerning 
the relationships with the quantum theory of angular momentum and group 
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representation properties of 50(3). We start from the analytic expression 
(valid for mi + m 2 + m 3 = 0, see [23], expression (8.2.1.5)) 



^1 h h 
mi ^2 m 3 



i)' i+mi V2^+T 



(h + h-hY-ih-h + hWi-h + hy. 



1/2 



(^3 + m 3 )\(l 3 - m 3 )\ 



ih + h+h + iy. 

1/2 



(h + mi)!(Zi - mi)!(Z 2 + m 2 )\{l 2 - m 2 )\ 

{-l) z {h + k + mi- z)\{h -mi + z)\ 



z\{l 2 +h~h- z)\(l 3 + m 3 - z)\(h -l 2 -m 3 + z)\ ' 



where the summation runs over all z's such that the factorials are non- 
negative. This expression becomes much neater for mi = m 2 = m 3 = 0, 
where we have 

^1 h h 




(34) 



, for h+l 2 + h odd 



■1) 



'-l+'2-'3 



Kh+h+h)/2]\ r (h 

[(h+h-h)/m(h-i 2 +h)/m(-h+i2+k)/2\\ \ 

for li + l 2 + I3 even 



+h-l3)Kh-h+l3)K-h+h+l 3 V- 



1/2 



Some of the properties to follow become neater when expressed in terms of 
the so-called Clebsch-Gordan coefficients, which are defined by the identities 
(see [23], Chapter 8) 

^1 h h 



(35) 



mi m 2 —m 3 



(36) C\l™l hm2 = {-l) h - l2+m3 y/2k + l ( ^ 



^1 h h 
mi m 2 —m 3 

We have the following orthonormality conditions: 
/oy\ \ r film f~il'vn! 



°l °m > 



mi ,ni2 



(38) 



Lm 



Now recall the general formula ([23]) eqs. 5.6.2.12-13, or [13J, eqs 3.64 and 
6.46) 

Y hmi (x)...Yi nmn (x)dx 



s2 



4?r 



2L + 



T S Yl [ C hn 

Li...L n _3 Mi...M n _3 



Mi S-1L2M2 
m\l2m2 LiMilzm-i'" L n -3M n -3l n -irn n -i 



(39) x 



(4 7r )n-l 



✓~Oj1U (~tL2§ f-lln.0 

o «i0/20 L/ Li0Z30-" o L„- 3 0Z I , 



-4 



Two important special cases are provided by 

*Jmi (^)^Zm 2 (a;)^m 3 fa)^ 
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_ , 1 \l-m 3 „ I fit + 1) n lQ n l-m 3 
' ' V 4^ ^I0l0^lmilm 2 



(40) 

and 



(21 + 1) 3 / I I I \ ( I I I 



4ir \ / V mi m2 m% 



f 



Y lmi {x)Y lm2 (x)Y lm3 {x)Y lm4 {x)dx 



LM n L ~ M 



(2l + l)y-^ \L—M (riLO l 2 Qmiim 2 ^im 3 im4 

] 1 im) 2TT1 



(41) 



I I L\ 2 ( I I L \ ( L I I 



L 

Similarly, the following identities hold: 

III 




I V mi ?7i2 M J \ M m% 7714 



L 



1 / 7 7 7 \ 2 

3/ 



/ I I L x 2 



Pf(t)dt 

o 

,i 2 * 
/ J?(t)dt = £(2£ + l) 

■ / ° L=0 V 

Finally, we recall the useful identity, valid for all x±,X2 £ S 

47T ^ — 

(42) P/((x!,x 2 )) = ^ Y lm ( Xl )Y lm (x2), 

m=—l 

which allows to express Legendre polynomials in terms of spherical harmon- 
ics. 





L=0 

2 
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